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Abstract 

We study hydrodynamic fluctuations in a non-relativistic fluid. We show that in three dimensions 
fluctuations lead to a minimum in the shear viscosity to density ratio rj/n as a function of the 
temperature. The minimum provides a bound on n/n which is independent of the conjectured 
bound in string theory, n/s > ^/(47tA;b), where s is the entropy density. For the dilute Fermi 
gas at unitarity we find n/n > 0.3h. This bound is not rigorous, but it can only be avoided if the 
domain of validity of hydrodynamics is anomalously small. We also find that the viscous relaxation 



time of a hydrodynamic mode with frequency u diverges as 1/y/uJ, and that the shear viscosity in 
two dimensions diverges as log(l/cj). 
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I. INTRODUCTION 



It is now widely appreciated that fluid dynamics can be viewed as an effective long- 
distance theory for a classical or quantum many-body system at non-zero temperature. 
Effective theories make systematic predictions for correlation functions order-by-order in 
a low-momentum expansion. These predictions depend on a small number of microscopic 
parameters. In the case of fluid dynamics the microscopic parameters are the equation of 
state and the transport coefficients. 

Effective (field) theories are constructed using the following procedure: i) Identify the low 
energy degrees of freedoms, ii) Write down the most general local effective action consistent 
with the symmetries of the problem. This action is typically expressed in terms of low 
energy fields and their derivatives. The coefficients of allowed terms in the effective action 
are free parameters called low energy constants, iii) Determine what terms in the effective 
action have to be included in order to compute a correlation function to a given order in 
the low energy expansion. This is known as the "power counting". Typically, the leading 
contribution arises from tree level diagrams involving operators with the minimal number 
of derivatives, and higher order corrections arise both from higher derivative operators, and 
from loop diagrams with the leading operators. In some cases diagrams may have to be 
summed to all orders. For example, the sum of all tree diagrams corresponds to solving a 
classical field equation. 

In fluid dynamics the low energy modes are fluctuations of the conserved charges. In the 
case of a one-component non-relativistic fluid the conserved charges are the particle density, 
the energy density, and the momentum density. The derivative expansion is implemented 
at the level of the constitutive equations. This means that the conserved currents are 
expressed in terms of derivatives of the thermodynamic variables. The effective theory 
without derivative terms in the currents is called ideal fluid dynamics, and the equation 
of motion at one-derivative order is known as the Navier-Stokes equation. The validity of 
fluid dynamics requires that derivative corrections to the currents are small. This condition 
does not preclude the possibility that small corrections can exponentiate as one solves the 
equations of motion. Solutions of the Navier-Stokes equation are qualitatively different from 
solutions of ideal fluid dynamics. In ideal fluid dynamics the motion is time reversible, sound 
modes are not damped, and diffusive modes do not exist. This implies that in most cases 
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one has to retain at least one- derivative terms in the constitutive equations. Two derivative 
terms have also been studied , but the corrections are typically small. In relativistic 
fluid dynamics two-derivative terms improve the stability of equations of motion, and second 
order terms are now routinely included in hydrodynamic simulations of relativistic heavy 
ion collisions [3j. 



uctuations are 
4(, but they are 



Loop corrections in fluid dynamics arise from thermal fluctuations, 
known to be important in the vicinity of second order phase transitions 
rarely considered in the case of non-critical fluids. In this work we will consider the contribu- 
tion of fluctuations to the correlation function of the stress tensor in a simple non-relativistic 
fluid. We study the implications of our results for the shear viscosity of the unitary Fermi 
gas. The unitary Fermi is known to have a very low viscosity js - 8 ] , close to the value ob- 
tained from the AdS/CFT (Anti-deSitter Space/Conformal Field Theory) correspondence, 
r]/s = h/(4irkB) @, Q- Here, r\ is the shear viscosity and s is the entropy density, % is 
Planck's constant and ks is Boltzmann's constant. We will set H = ks = 1 in the following. 
We will show that the small shear viscosity enhances the role of fluctuations, but we also 
show that fluctuations imply a lower limit on how small the viscosity can get. We demon- 
strate that fluctuations lead to a non-analytic term in the viscous relaxation time in three 
spatial dimensions, and to a logarithmic divergence of the shear viscosity in two dimensions. 
Finally, we discuss the possibility of observing these non-analytic terms in experiments with 
trapped atomic gases. 

Our work builds on a substantial literature related to fluctuations in fluid dynamics, 
beginning with the work of Landau [ll]. The role of fluctuations in critical transport phe- 
nomena was summarized in the review article by Hohenberg and Halperin |4J] and the text 



relativistic fluids, see 



17 



books by Ma and Onuki 
authors, for example in 14 



131 ] . Diagrammatic methods are discussed by a number of 
y follows a recent study of fluctuations in 
191 ] . and the related work in 20]. 



161 ] . Our work close 



18], the recent review 



II. KUBO FORMULA 



In this section we will determine the low energy behavior of the retarded correlation 
function of the stress tensor using the equation of fluid dynamics at next-to-leading order in 
the gradient expansion. This result can be used to derive the standard Kubo formula for the 
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shear viscosity, as well as a new Kubo formula for the viscous relaxation time. We will employ 



the formalism developed in 



QBE 



23], which is based on coupling the theory to a non-trivial 
background metric gij(t,x). Correlation functions of the stress tensor can be computed by 
using linear response theory, and the requirements of Gallilean and conformal symmetry can 
be incorporated by requiring the equations of fluid dynamics to satisfy diffeomorphism and 
conformal invariance. 

The retarded correlation function of the stress tensor IP- 7 is defined by 

G% kl (u,k) = -% J dt J rfxe^- ik - x e(t)([n^'(t,x),n fc '(o,o)]} . (i) 

Gr determines the stresses induced by a small perturbation gij(t, x) = 5ij + hij(t, x) around 
the flat metric. We have 

5W = - l -Gf l h kl . (2) 

In fluid dynamics we expand the stress tensor in derivatives of the local thermodynamic 
variables P, p, v, where P is the pressure, p is the density, and v is the fluid velocity. We 
write Uij = II - + SUij, where 

11°. = pv l v j + Pg tj (3) 

is the ideal fluid part, and 5Uij is the viscous correction. In a conformally invariant fluid 
the leading term is 8Hij = — rjaij with 

2 

<y%j = ViVj + VjVi + gij - -gij(cr) , (4) 
(a) = V-v + f g , (5) 

where is the shear stress tensor, 7] is the shear viscosity, gij(cr) is the bulk stress tensor, 
and Vj is the covariant derivative associated with g^. Note that the bulk viscosity of a 



conformal fluid is zero. In \l\ we classified all terms up to second order in derivatives. We 
have 

5^ = -rjaij + 7]t r (g ik a k j + v k V k Vij + ^(a)^ + ^ {l k ^ j)k + ^( k ^j) k 

+ \zn { k n j)k + ^v^tv^t + 72 v (4 pv i> p + l3 v { /rv j} p 

+ 74V (i V j} T + 7 5 V (i V,)P + K R R {ij) , (6) 

where tr is the viscous relaxation time, Ai, ji, kr are second order transport coefficients, 
&ij = VjWj — VjVi is the vorticity tensor, T is the temperature, and Rij is the Ricci tensor 



associated with gij. Note that Rij vanishes in flat space gij = Sij, but keeping terms involving 
the curvature is crucial for obtaining the correct low energy expansion of Gr. 

We will concentrate on the "pure shear" component G R yxy . For this purpose we consider 
a perturbation of the form h xy (z, t). From the linearized Euler equation we can see that the 
perturbation does not induce a shift in the density, temperature, or velocity. This means 
that we can directly read off 5Uij from equ. (Q. We find 

G x R y3:y {cj, k) = P - irjuj + T R r]co 2 - ^k 2 + 0(u 3 , uk 2 ) , (7) 

which implies the familiar Kubo relation for the shear viscosity 

77 = — lim lim 4— im G^ xy (uj, k) (8) 

as well as a new Kubo formula for the viscous relaxation time 

1 d 2 

trT] = lim lim -— Re G x R yxy {u, k) . (9) 
ui^O k^o 2 du z 



This result is simpler than the corresponding formula in relativistic hydrodynamics 
which also involves a term proportional to k,r. 



III. HYDRODYNAMIC FLUCTUATIONS 



In this section we will study the contribution of fluctuations to the retarded correlation 
function. For this purpose it is convenient to start from the symmetrized correlation function 

G x s yxy (u,k) = Jd 3 xJ dte^-*^ (±{n xy (t,x),n xy (0,0)}} . (10) 

This function is related to the retarded correlator by the fluctuation dissipation theorem. 
For u — > we have 

2T 

G s (u,k)~ ImG fl (w,k). (11) 

CO 

In the low frequency, low momentum limit we can use the form of the stress tensor in fluid 
dynamics, 11^ = pv x v y — rj(V x v y + V y v x ) + 0(V 2 ), and expand the hydrodynamic variables 
around their mean values, p = po + 5p etc. We will use the Gaussian approximation and 
write expectation values of products of fluctuating fields as products of two point functions. 
The ideal (zero derivative) terms in the stress tensor give one and two loop graphs involving 
velocity-velocity and density-density correlation functions. We will show in the appendix 
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that additional loops are suppressed by powers of 1/ (ra 1//3 / m ^,), where n is the number density, 
and l m f p is the mean free path. We will therefore concentrate on the one-loop graph 



r du' 


f d 3 k r 


1 27 J 


(2tt) 3 i 



G7>,o) = pS 



where A^ is the symmetrized velocity correlation function 



Af(u\k)Af(u-u\k)+Af(u\k)Af(u-u\]t)\ . (12) 



A?(w,k) = ^e^- kx )^p(t,x),^(0,0)}) . 



(13) 



We are ultimately interested in the retarded, not the symmetrized, correlation function. At 
low frequency the retarded function can be written as 



G2">,o) = pI 



r dJ 


r d 3 k r 


1 2^J 


(2tt) 3 i 



AgV, k)A^(w - w', k) + Af(oj f , k)A^(w - a/, k) 



+ Af(u/, k)A^(w - u/, k) + AfV, k)A^(w - to', k) 



(14) 



where we have used the fluctuation dissipation relation ( TTTT) . This relation generalizes: 
retarded correlation functions of hydrodynamic variables have diagrammatic expansions in 
terms of retarded and symmetrized correlation functions 
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19] . 



The velocity correlation function can be decomposed into longitudinal and transverse 
parts 

(15) 



A% R (co, k) = - W) Al R (u, k) + WA L StR (u;, k) 



The transverse part is purely diffusive. The symmetrized correlation function is 11] 

2T n u2 



A>,k) 



\2 ' 



(16) 



p io 2 + (D v k 2 y 

where k — \k\ and D v = rj/p is the momentum diffusion constant, also known as the kinetic 
viscosity. The retarded correlation function is given by 

„2 



A£(w,k) 



1 —D v k 



(17) 



p —iuj + D v k 2 

The longitudinal correlation function can be reconstructed from the density-density corre- 
lation function (the dynamic structure factor) using current conservation, iu5p = ipk ■ 
where v = vi + vy with k • vy = and k x \ L = 0. We find 



A L s (u;,k) 



2T 



Tuk 2 



+ 



P {(io 2 -c 2 k 2 Y + (Tiok 2 Y \c v 



Dtuj' 1 



c 2 uj 2 + (D T k 2 y 



- 1 

c v 



(tu 2 - c 2 s k 2 )D T uj 



c 2 s (u) 2 - c 2 k 2 y + {ruk 2 y 



(18) 



The first two terms have a clear physical interpretation as the contributions from propagating 
sound waves and diffusive heat transport. The third term is required to satisfy sum rules. 
This term is suppressed near the sound pole u 2 ~ c 2 s k 2 . In equ. (ITS]) c s is the speed of sound, 
T is the sound attenuation constant, D? = n/(c p p) is the thermal diffusion constant, k is 
the thermal conductivity, c p is the specific heat per unit mass at constant pressure and c v is 
the specific heat at constant volume. The sound attenuation constant is 
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3p p 
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= 3p _ 



3C , 3Ac p 



1 + -- + 



(19) 



4 77 4 Pr 

where Ac p = (c p — c v )/c v and Pr = (c p t))/k is the Prandtl number, the ratio of momentum 
to thermal diffusion. At high temperature Ac p = 2/3 and Pr = 2/3 25), and at low 
temperature Ac p /Pr — > 0. In the case of a conformal fluid the bulk viscosity vanishes, 
( = 0. This implies that T = \D r] at high temperature, and T = ^D v at low temperature. 

At low frequency and momentum the symmetrized correlation function can be further 
simplified. We illustrate the result in the case of the sound pole. We can write 



VTh 2 1 1 

AT nd (u,, k) * iif —r, + 



2 P {(u-c s k) 2 + (^y (co + c s kf + (!£)' 
which is correct up to terms of order Tk/c s . The retarded correlator is 



(20) 



^V.k)^) — ^ + — i^ }, (21) 

and an analogous expression holds for the sum rule term in equ. ([TBI) . 

We can insert the decomposition of the velocity correlation function given in equ. ( lisp 
into the one-loop result for the retarded correlation function, equ. (fl4|) . This gives a number 
of terms, corresponding to the contribution of shear modes, sound modes, a mixed shear- 
sound contribution, and the contribution of diffusive heat modes, see Fig. [TJ We discuss 
these contributions in turn: 

1. Shear modes: This is the contribution which is easiest to compute. The frequency 
integral can be done by contour integration. We find 

7T r k 4 

This integral is UV divergent. We regulate the divergence by introducing a momentum 

cutoff A. We then expand the integral in the low frequency regime. We get 

7 7TA IT 

<*■>• 0)1- - - jjjj^A' - iu, ^ + (1 + + O(^) . (23) 



a) ^ ^ b) r y^^ c) 




o- 



FIG. 1: Diagrammatic representation of the leading contribution of thermal fluctuations to the 
stress tensor correlation function. Solid lines labeled vt denote the transverse velocity correlator, 
dominated by the shear pole, and wavy lines labeled vl denote the longitudinal velocity correlator, 
governed by the sound pole and the diffusive heat mode. 

The physical meaning of these terms can be understood by comparing with the Kubo relation 
in equ. (JZJ). The first term is a fluctuation contribution to the pressure, and the second 
term is a correction to the shear viscosity. The imaginary part of the third term can be 
viewed as a frequency dependent correction to the shear viscosity, and the real part is a 
frequency dependent contribution to the relaxation time which diverges as uj~ 1 ^ 2 in the low 
frequency limit. The existence of this term is sometimes interpreted as an indication that 
hydrodynamics breaks down beyond the Navier-Stokes (one-derivative) order. 

2. Sound modes: In order to calculate the contribution from sound modes we use equ. ( 12 Op 
and ( l2Tj) . This leads to two types of terms, depending on whether the real parts of the poles 
of the propagators in the w-plane have the same or opposite sign. The contribution from 
terms with opposite real parts has the same structure as the shear mode term. We get 

o)L- = ^ - * + a + i)^^r 2 + o ( ^) . (24) 

The contribution from terms with real parts of the same sign is not infrared sensitive and 
does not contribute to the retarded correlation function at O(co) or 0(tu 3 ^ 2 ). 

3. Shear-sound contribution: The shear-sound contribution has the structure 

CT.°)I,W> p + i( ^m - (25) 

where D s = D v + T/2 and the range of the /c-integral is [—A, A] because of the k -H- —k 
symmetry of the sound propagator in equ. fl20l) . We get 

GT V (^ 0)L_ ~ ^ + + 0(c 2 ) , (26) 



which is suppressed relative to the pure shear and sound contributions by a factor D s A/c s <C 
1 (see Sec. IIVI) . We also note that the mixed shear-sound term does not give non-analytic 
terms of the form w 3//2 . 

4. Diffusive heat modes: The contribution of diffusive heat modes is very similar to the 
shear term, but the residue of the heat mode is proportional to u 2 /c 2 instead of k 2 . In the 
diffusive regime u 2 <C c 2 k 2 . We find 

Gr> , 0)U „ !^/*_*_ . <**L, (27) 

which is much smaller than the shear term. 

We conclude that the main contribution arises from the pure shear and sound terms. We 
will combine these two contributions using the approximation T ~ ^D v , which corresponds 
to the low temperature regime. This is the more interesting regime because D v is small and 
the role of fluctuations is enhanced. We find 

17TA 7+(~Y 
G2~( U , 0)L = const - iu, + (1 + ^T^Lj. + o^) . (28) 

As noted above the iu term is a contribution to the shear viscosity. This term is cutoff 
dependent, but the physical viscosity must be independent of an arbitrary cutoff. This 
implies that the bare viscosity must be cutoff dependent too, and that the cutoff dependence 
of the bare viscosity is governed by a renormalization group equation. It is important that 
the non-analytic co 3 ^ 2 term is not cutoff dependent, because any cutoff dependence in this 
contribution cannot be absorbed into the parameters of the hydrodynamic description. 

IV. PHENOMENOLOGICAL ESTIMATES 

In this section we study phenomenological implications of the results derived in the previ- 
ous section. We observed that iu term in the retarded correlation function can be combined 
with the bare shear viscosity to give a physical viscosity 



17 pTA 



An interesting consequence of this result is the fact that the physical viscosity cannot be 



arbitrary small 17|, because equ. ( 129|) has a minimum as long as the bare viscosity is positive. 



The bare viscosity must be positive for the hydrodynamic expansion to be well defined. The 



value at the minimum depends on the value of the cutoff; the larger the cutoff the stronger 
the bound on rj becomes. The largest possible value of the cutoff is determined by the 
condition that the gradient expansion on which hydrodynamics is based must be valid for 
all k < A. In the following we will study this condition separately in the shear and sound 
channel. 

1. Shear channel: Shear modes are characterized by u ~ D v k 2 . Corrections arise from 
higher order terms in the derivative expansion. For non-zero frequency the leading correction 
is due to the relaxation time. We have tu ~ D v k 2 <C . For this relation to be maintained 
for all k < A we need to require that A < K hyd with K hyd = (trD^ 1 / 2 . In kinetic theory 
tr = v/P Q, 23, 26] and 



1 fP\ 1/2 

^^(-j . (30) 
2. Sound channel: In the sound channel we have lo ~ c s k <C Tk 2 . Using T ~ we find 

^i©.* (31) 

For a weakly interacting gas (dP) / (dp) s ~ (5P)/(3p), and equ. fl30|) differs from equ. fl3Tj) 



by a factor very close to one, y 16/15 ~ 1.03. In the following we will use equ. ( 1301 as our 
estimate for the cutoff. We note that near a critical point the speed of sound can go to zero, 
and the contribution of sound waves is strongly suppressed relative to shear modes. 

It is interesting to consider the microscopic meaning of the ultraviolet scale Kh y d- In 



kinetic theory t] ~ npl m f p , where p ~ y/rriT is the mean momentum, and l m f p is the 
mean free path. For a weakly interacting gas P ~ nT and the ultraviolet scale is 
Khyd = {p/v){P/ P) 1 ^ 2 ~ Cjp- This is physically reasonable: It does not make sense to 
consider hydrodynamic fluctuations with wavelengths shorter than the mean free path. 

Finally, we can estimate the minimum shear viscosity. Using equ. ( )29|) and A = 
from equ. fl30|) we get 

, - ■ , , , ^ 

,nj phys n ^7280 \T F J \r] 



") = ? + JLm ! "W f^" 2 (32) 
a/2 80 \T f J \v \nTJ ' v ' 



where we have introduced the Fermi temperature Tp = k 2 F j (2m) and the Fermi momentum 
k F = {S^n) 1 ^. Minimizing this expression with respect to the bare viscosity we find 

3 /289\ 1/6 / T \ 1/2 / P \ 1/6 / T \ 1/2 ( P \ 1/6 , . 
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(-) 



T/T 



T/T 



FIG. 2: Shear viscosity to density ratio rj/n as a function of T/Tp, where Tp is the local Fermi 
temperature. The left panel shows r\/n for the unitary gas in three dimensions. The solid line is the 
result in kinetic theory and the dashed line includes fluctuations. The band shows the uncertainty 
if the cutoff is varied in the regime A = (0.5 — 1.0) Kh y d . The right panel shows the two-dimensional 
gas at the crossover point T a ^d = Tp. The solid line is the kinetic result. The dashed and dotted 
lines include fluctuations where we have used A = K^yd and u) was taken to be the frequency of 
the quadrupole mode in a harmonic trap with N = 10 4 and N = 10 5 particles. 

Fluctuations imply a bound on the shear viscosity that decreases with the temperature 1 . 
Kinetic theory calculations show that the mean free path and the bare viscosity become 



very large in the limit T <C T c [271, |28[. It is therefore most interesting to consider the 
implications of the bound near the critical temperature T c . In this regime 
Using T c = 0.167(13)Tp (29 1 we get (r]/n)T=T c ?Z 0.41. There is some uncertainty is this 
result due to higher order loop corrections, and due to the uncertainty in the cutoff. We note, 
however, that the bound scales as A 1 / 3 , so that a factor 2 uncertainty in A only translates 
into a 25% uncertainty in the bound. As a conservative estimate we take A ~ O.bK^ and 

V 



> 0.32. (34) 

n J t=Tc 

We can illustrate this result further by using the leading order kinetic theory result as an 
estimate for the bare viscosity. This is consistent because kinetic theory takes into account 



1 We note that the lower bound on rj/s is attained at a non-zero temperature. In particular, (rj/s) p h y s > 
1.005(n/s)(T/T F ) 1 / 2 (P/(nT)) 1 / 6 is minimal at a temperature T - T c becausejs/n) - T 3 for T < T c 
and (s/n) ~ log(T) for T ^> T c . Using the entropy per particle measured in [29( and A = 0.5Kh y d, we 
obtain a fairly broad minimum (rj/ 's) p hy S > 0.3 in the regime T /Tp ~ (0.3 — 0.8). 
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effects at distances / < l m f p but, unless additional resummation is performed, it does not 
include fluctuations at length scales I > l m f p . The kinetic result is [30] 

^^T?n. (35, 

In Fig. [2] we show the bare viscosity and the physical viscosity including the effects of 
fluctuations. We observe that the viscosity has a minimum rj/n ~ 0.5 at a temperature 
T ~ 0.27>. 

At the phase transition the entropy per particle is s/n ~ 0.8 and equ. f[3"4"l) translates into 
{t)/s)t c <; 0.4. This number is comparable to experimental results reported by Cao et al. 6| 
and to the T-matrix calculation of Enss et al. jsij], but larger than the Path Integral Monte 
Carlo (PIMC) results obtained by Wlazlowski et al. [32|. These authors find (i]/s) min > 0.2, 
quite close to the AdS/CFT bound (rj/s) > l/(47r). If the low value of rj/s obtained in the 
PIMC calculation is indeed correct then the breakdown scale of hydrodynamics would have 
to be much smaller than the estimate given in equ. ( 130|) . Another possibility is that for 



the lattice spacing used in 



32] one cannot resolve non-analytic contributions to the spectral 



function. Equ. ( 128]) implies that for small u 



3/2 



7+ faV 

V (u) = V -V^T ^3^. (36) 

This relation shows that for small values of r\ it is difficult to resolve the zero frequency limit 
of the shear viscosity. The width of the non-analytic structure in the spectral function can 
be estimated by assuming that the fluctuation term in equ. (|29|) is related to the y/oJ term 
in the spectral function. This assumption leads to Au ~ 0.3T(n/?7), where we have used 
A ~ 0.5K hyd . 



V. TWO DIMENSIONAL SYSTEMS 



It is interesting to consider the role of fluctuations in two dimensional systems. The 
results in Sect. IIHI are easily generalized to two spatial dimensions. Aside from the obvious 
substitution d 3 k/(27r) 3 — > d 2 k/(2ir) 2 the only difference is that in two spatial dimensions 
the shear contribution to the sound attenuation constant is T = D v instead of T = ^D v . In 
both two and three dimensions the dominant contribution to G^ xy arises from the one loop 
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diagram involving either a pair of shear modes or a pair of sound modes. In d = 2 the loop 
integral is logarithmically divergent and 



I tot 



const — iu ■ 



T 



16nD r 



log 



{ V2D V A 2 \ ,7T 



(37) 



The imaginary part can be interpreted as a correction to the shear viscosity. We find 

mT n ( \^2A 2 r\ \ 



167T T) 



log 



\ mw n 



(38) 



which diverges logarithmically as u — > 0. This divergence is well known 33M35| . and it has 



been observed in molecular dynamics and lattice gas simulations [36|, [37J . To the best of our 
knowledge it has not been observed experimentally. Equ. (1381) shows that the effect is large 
in systems that have a small value of the bare shear viscosity. 

The shear viscosity of a dilute two-dimensional Fermi gas was recently computed in 

log 



38l . |39j . The result is 



mT 



V 



2vr 2 



5T 



n 2 



(39) 



where T a ^d = l/(^ia2d) 2 and a,2d is the scattering length in two dimensions 40]. In two 
dimensions there is no scale invariant fluid except in the non-interacting limit a 2 d — > 0. 
The most strongly correlated fluid corresponds to T a ^d — Tp, which implies that the dimer 
binding energy is equal to the Fermi energy. The viscosity of the two dimensional gas was 
recently studied by measuring the damping of the quadrupole mode [4l| in a harmonic 
potential. The frequency of this mode is u = y/2u±, where u± is the two dimensional 
oscillator frequency. The confinement frequency sets the scale for the Fermi temperature of 



the trap, T F 



trap 



We can use these relations to translate the frequency dependence 



of the shear viscosity into the dependence on the number of particles. In Fig. [2] we show the 
kinetic theory result as well as the viscosity with fluctuations included for two different values 
of the particle number. The main dependence on N is of the form (r)/n) p h ys ~ ^ log(iV). 
We observe that fluctuations make a significant contribution to the shear viscosity, but the 
logarithmic divergence with iV is fairly slow, and one will need significantly larger numbers 
of particles than what is available in current experiments (N ~ 4 • 10 3 in [411]) to see the 
effect clearly. 



13 



0.0 0.2 0.4 0.6 0.8 1.0 

T/T 



FIG. 3: Trap averaged shear viscosity to density ratio (a n ). We show (a n ) as a function of T/Tp ap , 
where Tp ap = {"iXNY^ui^ is the Fermi temperature of the trap. We have chosen N = 2 ■ 10 5 and 
A = 0.045 as in [43J. The solid line shows the kinetic theory result, the dashed line includes 
fluctuation corrections to the shear viscosity. The data are from [44 |. which is a reanalysis of the 



results reported in 



43J. We do not show data in the superfluid regime T -C T c . 



VI. TRAPPED ATOMIC GASES 

In this section we will try to make contact with experiments that study the damping of 
collective modes in trapped Fermi gases. We are interested in the question whether it is 
possible to establish the role of hydrodynamic fluctuations by studying the scaling of the 
damping constant with temperature or particle number. A review of the hydrodynamic 



theory of collective modes can be found in 



42] 



Consider a trapped gas with N particles in a harmonic potential with trapping frequencies 
uu = uj y = cj± and u z = \u±. In a typical experiment N = (10 5 — 10 6 ) and A = (0.02 — 0.05) 



43 



45 



-|47j. The transverse breathing mode has a frequency u = Jl0/3co± and the damping 



constant is 

{a n ) w± (m) 

Lbr {W\y/*{Eo/[Ne F ])' 1 } 

where Eq is the total (potential and internal) energy of the trapped gas, = (SNXy^tu^ 

is the Fermi energy of the trapped system, and (a n ) = f d 2 xrj(x) is the trap average of 
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the shear viscosity. Taking into account relaxation time effects we have 



(«n) = 7jffc 1 7 7 ,2 • ( 41 ) 

A' J 1 + u z t r {xY 



We take the bare shear viscosity from kinetic theory, equ. (|35p . and compute the physical 
viscosity from equ. fl32|) . We note that the bare viscosity only depends on T, which is 
independent of the position in the trap. The fluctuation term is largest at the center of the 
trap. We also use kinetic theory to determine the bare relaxation time, tr = r]/P, and use 
equ. ([25]) to determine the physical relaxation time. This corresponds to 

3/2" 



7 +(D 

(TRV) P kys = T RV + 24(W/2 " 



T ( 



(42) 



WJ 

We note that the bare relaxation time is inversely proportional to the local pressure and 
depends on the position in the trap. In particular, trT] is large in the dilute part of the 
cloud. Fluctuations, on the other hand, increase the relaxation time near the center of the 
trap. 

In the following we will use the high temperature approximation for the density of the 
cloud. This is consistent with using kinetic theory for the bare shear viscosity and relaxation 
time. It also provides a very accurate description of the tail of the density distribution at 
essentially all temperatures. We have 

"(^gfH?^)' (43) 

where u> = (ur[u; 2 ) 1//3 . Results for (a n ) as a function of T/T F rap with Tp ap = (SiVA) 1 / 3 ^ are 
shown in Fig. El We have used N = 2 ■ 10 5 and A = 0.045 as in the experiment of Kinast 



et al. 43(. The solid line shows the result using kinetic theory for rj and tr, and the dashed 
line includes the fluctuation term in r]. We find that for the parameters considered here 
corrections to the relaxation time are very small. 

We observe that kinetic theory describes the data for T > 0ATp ap well. Fluctuations 
are important for T < 0.2Tjr ap , leading to a minimum in (a n ). We note that the critical 
temperature is T c ~ 0.2Tp ap [48 1. and we do not expect the theory used in this section, 
which is based on kinetic theory in the dilute limit, to reproduce experiment for T T c . 
It was recently suggested that the data in this regime are dominated by the transition from 



hydrodynamic to ballistic behavior 



42| 
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Finally, we have looked at the role of fluctuations in the experiment of Vogt et al. [41 ]. 
In this experiment the damping of two dimensional quadrupole mode was measured for 
N = 4 ■ 10 3 . The dependence on the scattering length was studied for \og{kpa 2 d) > at 
T/Tp = 0.48, and the temperature dependence was studied in the range T/T F rap = (0.3—0.8) 
for \og(kF(i2d) = (2.7 — 42). We find that for this range of parameters the role of fluctuations 
is always small. Fluctuations lead to a significant enhancement of the damping constant if 
log(kpa 2 d) ~ and iV > 10 5 . This enhancement grow as log(iV), but the logarithmic growth 



is difficult to disentangle from a log(iV) term related to the dilute corona, see [38 1. 



VII. CONCLUSIONS AND OUTLOOK 

We have studied the role of hydrodynamics fluctuations in the dilute Fermi gas. Our 
main findings are: 

1. Hydrodynamic fluctuations lead to a minimum in the shear viscosity. The physical 
origin of the minimum is the contribution of shear and sound modes to momentum 
transport. If the bare viscosity is small sound and shear modes are weakly damped 
and this contribution becomes large. The magnitude of the minimum shear viscosity 
is weakly sensitive to the cutoff scale for the hydrodynamic description. Allowing for 
a factor of two uncertainty in our estimate of A we find rj/n > 0.3. This bound is not 
rigorous, but it can only be violated if the domain of applicability of hydrodynamics is 
anomalously small. Our estimate is close to the trap averaged value of rj/n at T ~ T c 
measured experimentally p , but larger than the minimum rj/n ~ 0.15 obtained in a 
recent lattice simulation [32] . We note that the bound itself is purely classical, fi only 
enters though thermodynamic quantities 2 . 

2. Fluctuations cause a 1/ \fui divergence of the viscous relaxation time in three dimen- 
sions, and a log(u;) divergence of the shear viscosity in two dimensions. These effects 
are independent of the cutoff and only depend on the value of the bare shear viscosity. 



2 Consider the bound on -q/s. Ignoring numerical factors we have rj/s > (n/s)(mT/n 2 / 3 ) 1 / 2 (P/(nT)) 1 / 6 . 
At large temperature the ratio n/s depends weakly on T and the bound grows as T 1 / 2 . At low T the 
entropy per particle increases sharply when the system reaches quantum degeneracy, which corresponds 
to mT ~ h 2 n 2 / 3 . This implies that at the minimum rj/s ~ h. 
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The existence of non- analytic terms implies that, strictly speaking, the two dimen- 
sional Navier-Stokes equation as well as the three dimensional second order (Burnett) 
equations are not consistent unless fluctuating forces are taken into account. We note, 
however, that real flows that can be studied in experiment involve finite frequencies 
or time scales, and fluctuating forces may not be important. 

3. We have studied the importance of fluctuations for the damping of trapped Fermi 
gases. The corrections are generally small in the range of parameters that have been 
experimentally investigated. A possible exception is the three dimensional unitary gas 
near T c . In this case fluctuations may lead to enhanced damping 3 . Fluctuations lead 
to a log(iV) divergence in the damping constant of the two dimensional Fermi gas, but 
this effect is difficult to observe unless the number of particles is varied by more than 
an order of magnitude. 

There are a number of interesting formal questions that we have not studied in this paper. 
In order to study higher order corrections it is useful to start from an effective action for 
hydrodynamic fluctuations. This could be done using the methods developed in [lfl fioj ] . 
The effective action might also be useful for studying the renormalization group evolution 
in more detail. In this work we have simply assumed that the bare shear viscosity can be 
computed in kinetic theory. It would be desirable to provide a more rigorous justification for 
this approximation by studying the matching between kinetic theory and hydrodynamics. 
Finally, in this work we have restricted ourselves to studying the effect of fluctuations on 
the damping of collective modes. In this case it is straightforward to take into account a 
frequency dependent shear viscosity and relaxation time. If one considers hydrodynamic 
flows that are not periodic in time, for example the elliptic flow experiment described in jf|, 
one has to solve the hydrodynamic equations with fluctuating forces. This method has been 
studied in the context of microfluidic systems, see 50| and references therein. 

We have shown that hydrodynamic fluctuations are important if the bare viscosity is 
small in the low temperature limit. The main physical question is whether this scenario is 



3 Note that we have not considered the role of critical fluctuations. The superfluid transition is described by 
model F in the classification of Hohenberg and Halperin Q • This model does not contain direct couplings 
between the order parameter and the momentum density, and the calculation discussed in our work is not 
directly affected by critical fluctuations. 
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realized in the unitary Fermi gas near T c , or whether other effects, like pairing correlations, 



a pseudo-gap or phonons are more important 



31 



5lN53| . This question is probably difficult 



to address experimentally, but it can be studied by analyzing the spectral function of the 
stress tensor. For this purpose it is necessary to construct models of the spectral function 
that include fluctuations. These models can be confronted with quantum Monte Carlo data, 
for example the recent work of Wlazlowski et al. 32 ]. 
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Appendix A: Stress tensor correlation function in hydrodynamics 
1. Low energy expansion 

In this section we provide some additional details regarding the low energy expansion 
of the retarded stress tensor correlation function. Our starting point is the symmetrized 
correlation function 

G x s ' xy (u,k) = Jd 3 xJ dte^-^ (^{U xy (t,^),U xy (0,0)}^ . (Al) 

We use the expression for the stress tensor in hydrodynamics and expand in small fluctua- 
tions, 5p, 5T, Vi, and in the number of derivatives. We get 

n xy = p v x v y - i]o {V x v y + V y v x ) - 

where we have dropped terms of order 5 3 and V 2 . The diagrammatic expansion can be 
derived by inserting equ. (1A2|) into equ. (lAlj) and factorizing the expectation value into 
pairs of fluctuating fields. The retarded correlation function is obtained by replacing one of 
the symmetrized functions by a retarded function. In Section UTTl we computed the one loop 
diagram that arises from the first term in equ. ( 1A2I) . This term has no spatial derivatives, 
and we find a contribution of the form 

GT>, 0) ~ TA 3 1 1 + + + ..], (A3) 
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Sp + 




(V x Vy + VyV x ) + 



(A2) 




FIG. 4: Diagrammatic representation of higher order fluctuation contributions to the stress tensor 
correlation function. Solid lines labeled vt denote the transverse velocity correlator, wavy lines 
labeled vl denote the longitudinal velocity correlator. Temperature fluctuations are shown as 
double dashed lines, and density fluctuations are shown as dotted lines. Vertices shown as dots 
contain no derivatives, whereas vertices labeled by squares contain spatial derivatives. 

where C\ and C3/2 are numerical constants. We observe that the low energy expansion involves 
powers of u/^D^A 2 ). 

The second term contains spatial derivatives, and it gives rise to a tree diagram (Fig. H^) 
which vanishes as k — > 0. The third term in equ. (1A2I) involves derivatives of the shear 
viscosity with respect to p and T. In kinetic theory [(dr]) /(dp)]r vanishes at leading order 
in n\ 3 dB , where Xjb is the de Broglie wave length. The dominant term therefore involves 
fluctuations of the temperature, G$ ~ [(dr])/(dT)] 2 (5T5T){V x v y V x Vy) (plus permutations 
x — > y), see Fig.Hb. This is a one- loop graph with vertices that contain one spatial derivative. 
For non-zero external momenta this graph gives a contribution to Gr which is suppressed 
by k 2 /K 2 yd relative to equ. ( 1A3I) . For zero external momentum the diagram has power 
divergences which contribute at 0(A 2 / 'K 2 yd ). This is not small if A ~ Kh y d, but power 
divergences can be absorbed into the transport coefficients 4 . 

The leading no-derivative contribution of order 5 3 in the stress tensor is of the form 
n x j/ = (bp)v x v y . This term generates the two-loop diagram shown in Fig. Hfc. The extra loop 



It is well known that power counting is not manifest in effective field theories regularized by a momentum 
cutoff 55]. This problem can be circumvented using dimensional regularization (DR), which automatically 
eliminates all power divergent terms. In our context this is a disadvantage because we find that the leading 
one loop divergence in the stress tensor correlation function represents an important physical effect. The 
linear divergence is preserved in a modified version of dimensional regularization called "power divergence 
subtraction" (PDS) [scj , which keeps the pole corresponding to the logarithmic divergence in two spatial 
dimensions. 



19 



integral involves three powers of momentum, and the extra propagator is proportional to 
l/p. As a consequence, the graph is suppressed by A 3 /kp. This can be written as (A/fcp) 3 < 
(K h y d /kF) 3 — l/(n/^ p ), where we have used the kinetic theory estimate r] ~ n(mT) 1 / 2 / fn j p . 

In summary, the hydrodynamic expansion of the correlation function G°^(u,k) involves 
powers of k/K hyd , where K hyd = (P / p) 1 / 2 D~ l ~ is the breakdown scale defined in 
Section [TV] Since G^(u, k) is diffusive the frequency scales as u ~ D v k 2 and the frequency 
expansion involves (u/[D v K 2 yd \). Powers of k/Kh yd arise from higher derivative terms in 
the currents or from loop graphs. Loop graphs are additionally suppressed by (Khyd/kp) 3 — 
1 / {nl^fp) ■ Loop graphs are important because they lead to non-analytic effects, logarithms 
and fractional powers of u and k 2 . For low viscosity fluids the mean free path is short, 
n tfnh ~ 1> an d there is no suppression of loops relative to gradient terms. 



2. Contact terms 



The symmetrized correlation function contains a contact term which can be determined 



using the fluctuation-dissipation theorem 

1 



11 



57| 



2 {U xy (t,^),U xy (t',^')}j = 2r / T5(t-t')5(x-x / ) . (A4) 

In frequency space this gives the contribution of the bare shear viscosity to the retarded 
correlation function, Gr = iurj. The result therefore justifies combining the bare and 
loop corrections as in equ. (129]) . The contact term can also be obtained using the velocity 
correlation function combined with the conservation laws [17|]. Momentum conservation 
implies VjL™ = —-j^(pvi) and 

k 2 x GT y ^,k x )=co 2 Af(co,k x ). (A5) 
Using the explicit form of the velocity correlation function we find 

<??■>,*.) = ^{l-J^y}, (A6) 
which contains the contact term 2r]T. 
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